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C\ ! Abstract 

£h ■ We investigate the magneto-conductance (MC) anisotropy in the variable 

range hopping regime, caused by quantum interference effects in three dimen- 

o" _ 

O ■ sions. When no spin-orbit scattering is included, there is an increase in the 



localization length (as in two dimensions), producing a large positive MC. By 



contrast, with spin-orbit scattering present, there is no change in the local- 
ization length, and only a small increase in the overall tunneling amplitude. 
The numerical data for small magnetic fields B, and hopping lengths t, can 
be collapsed by using scaling variables B±t 3 ^ 2 , and Bnt in the perpendicular 
and parallel field orientations respectively. This is in agreement with the flux 
through a 'cigar'-shaped region with a diffusive transverse dimension propor- 
tional to y/i. If a single hop dominates the conductivity of the sample, this 
leads to a characteristic orientational 'finger print' for the MC anisotropy. 
However, we estimate that many hops contribute to conductivity of typical 
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samples, and thus averaging over critical hop orientations renders the bulk 
sample isotropic, as seen experimentally. Anisotropy appears for thin films, 
when the length of the hop is comparable to the thickness. The hops are then 
restricted to align with the sample plane, leading to different MC behaviors 
parallel and perpendicular to it, even after averaging over many hops. We 
predict the variations of such anisotropy with both the hop size and the mag- 
netic field strength. An orientational bias produced by strong electric fields 

will also lead to MC anisotropy. 
71.55.Jv, 05.40.+j, 72.20.Dp, 75.10.Nr 
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I. INTRODUCTION 



Striking quantum interference (QI) effects have been observed in experiments on insulat- 
ing materials |l]-f§|. These observations are of particular interest, as they point to quantum 
coherence phenomena for strongly localized electrons, where naively they may not have been 
expected to occur over length scales appreciably larger than the localization length £. A 
theoretical explanation was first proposed by Nguyen, Spivak and Shklovskii (NSS) || in 
the context of Mott Variable Range Hopping ( VRH) || : Phase coherence is maintained over 
the long distances between phonon assisted tunneling events, which grow with decreasing 
temperature T as exp(To/T) 1 /( jD+1 ^ in D spatial dimensions. The resulting coherence length 
can be quite large (typically of the order of 20 — 50£). In this work, we consider the three di- 
mensional NSS model and focus on the dependences of the conductance and its fluctuations 
on the relative orientations of the magnetic field and the dominant hop. 

The initial indications of QI came from observations of a strong positive magneto- 
conductance (MC) in materials that exhibit VRH behavior JF]. In a single impurity pic- 
ture, the action of the magnetic field is to further confine electrons already localized around 
the impurity. This would result in a negative MC which is not the case in experiments 
for weak magnetic fields. Further evidence is provided by the orbital nature of the MC 
observed in InO films of varying thickness. While experiments show an isotropic MC 
for thick samples, anisotropy sets in when the film thickness is close to the Mott hopping 
length. Such anisotropy precludes explanations in terms of scattering of electron spin by 
magnetic impurities |§, which are necessarily isotropic with respect to the field direction, 
pointing instead to interference effects due to the electron orbits. Finally, in a careful set 
of experiments, Orlov et al [||] and Milliken and Ovadyahu demonstrate the presence of 
reproducible conductance fluctuations or magneto-fingerprints, generally regarded as a clear 
signature of QI effects. 

The NSS model considers the QI between the many virtual paths that the electron can 
take while tunneling under the barrier between two distant impurity centers. In the tunnel- 
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ing process, the hopping electrons with energies near the Fermi level, encounter impurities 
with energies outside the Mott energy strip. These impurities are considered as the source 
of elastic scattering events under the barrier. Since the contribution of each virtual path 
(tunneling through a barrier) is exponentially damped by the distance it covers, it is suffi- 
cient to ignore back-scattering and focus on the (directed) paths that only undergo forward 
scattering between the initial and final impurities. 

The initial numerical studies of the NSS model (on relatively small systems) || indeed 
confirmed that it yields the correct sign for the MC. Subsequently, Sivan, Entin-Wohlman, 
and Imry || provided a theoretical analysis that agrees with much of the early NSS results. 
The critical hop is identified from the condition of producing a percolating network of random 
resistors |L0| , while the probability distribution for individual hops is calculated by assuming 
that the contributing virtual paths are uncorrelated. The latter assumption, which we 
shall refer to as the Independent Path Approximation (IPA), was shown to be invalid by 
Shapir and Wang ||12|| , since, in low dimensions, the paths must intersect at some scattering 



sites. Eventually, the correct form of the hopping probability distribution was calculated 
by Medina et al |13| , by incorporating the correlations between the virtual paths. The 



analytical results, confirmed by extensive numerical simulations on very large sizes, indicate 
that the positive MC in this model actually corresponds to an increase in the localization 
length with the magnetic field in the absence of Spin-Orbit scattering. This prediction is 
supported, at least qualitatively, by recent experiments on InO [[J and YSi [^4| . While the 



IPA scheme cannot produce a change in the localization length, an alternative approach to 
strong localization, based on Random Matrix Theory (RMT) JTJj], also produces such an 
effect. The latter, which is exact only for quasi-one dimensional systems, predicts a doubling 
of the localization length. 

There are conflicting theoretical and experimental observations in the presence of Spin- 
Orbit scattering. The first experimental study on InO [I5[] showed MC behavior resembling 
that of the weak localization regime; i.e. a positive MC for low fields, changing to negative 



at higher fields. On the other hand, more recent experiments |14| on YSi show negative 



MC for all applied fields. On the theoretical side, both the IPA scheme [ 16(1 , and the correct 
accounting of correlations ||17|| , yield a positive MC without changes in the localization 
length. (In fact the IPA results for MC are exact in the presence of strong SO.) Finally, 
the Random Matrix |14| approach finds a negative MC caused by a universal decrease of 
the localization length by a factor of two. These issues are discussed in greater detail in 



reference ||19|| , and will not be discussed further in the present work. 

In this work we study the NSS model in three dimensions, with and without spin-orbit 
scattering. An important new feature is that we must now take into account the relative 
directions of the magnetic field and the dominant hop. This issue is most relevant experi- 
mentally for samples that are small enough (or at such low temperatures) to include only a 
single dominant hop. By measuring the MC anisotropy as a function of the direction of the 
magnetic field, it is possible to locate the orientation of this dominant hop! Field depen- 
dences parallel and perpendicular to this orientation can then be used to further test the 
current models of coherence in the localized regime. There is, however, a certain amount of 
internal averaging when the conductance is dominated by several hops. Some insight about 
the nature of the hops can then be obtained by examining the MC anisotropy of thin films, 
as a function of their thickness and orientation to the magnetic field. 



II. THE NSS MODEL 

Low temperature conduction in the strongly localized regime is dominated by thermal 
hopping. At the lowest temperatures localized electrons lack enough thermal energy to hop 
to neighboring sites. On the other hand, electrons cannot wander too far away from their 
localization point due to the exponential decay of the wave function. The balance of these 
competing tendencies results in an optimal hopping length and leads to Mott's law for VRH 
||. Each of these hops may be represented by an effective resistor (hopping probability) 
in a network which can then be solved for the macroscopic conductance of the sample. 
The picture of the Miller-Abrahams (MA) |2D| network is central to the understanding 
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of hopping conduction. The effective resistance of the MA network was first estimated 



by Ambegaokar, Halperin, and Langer (AHL) |ID| , and Shklovskii and Efros ||11|| , using a 
percolation argument: Due to the exponentially large values of the resistors in the network, 
the macroscopic conductance is dominated by a single bottle-neck resistor on a percolating 
network. (We shall later discuss the modifications due to multiple hops). This simple 
argument provides a powerful tool, since it is then sufficient to determine the variations of 
a single hop with various external (applied fields) or internal (doping, correlation effects, 
anisotropy) physical parameters |2lJ. 

The model proposed by NSS examines QI effects for the dominating hop. Due to the 
long distance of the hop, typically R ~ ^(T /T) 1 ^ D+1 ) ~ 20 — 50£, electrons scatter off 
many impurities on route to the final site. While at the end of the process there is some loss 
of phase coherence (due to inelastic scattering by phonons), the intermediate scattering is 
elastic. To study QI processes for the hop, the NSS model places the impurities on the sites 
of a regular lattice; e.g. the cubic lattice in Fig. |T[ The interference effects are maximized if 
the initial and final sites for the hop are chosen at diagonally opposite end-points. Electrons 
can then follow many different virtual paths from the initial to the final site. The overall 
tunneling amplitude is computed by summing all (virtual) paths between the two points, 
each contributing an appropriate quantum mechanical complex weight. These weights are 
obtained from an Anderson tight-binding Hamiltonian 

H = ti<4 a i + H V ij a i a 3' (!) 

i <ij> 

where ei are the impurity site energies, and Vij represent the nearest neighbor couplings 
or transfer terms. NSS further simplify the problem by choosing site energies distributed 
according to 

{-\-W with probability p, 
—W with probability (1 — p), 

and a transfer term 

C V if i, j are nearest neighbors, 



otherwise. 
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We shall henceforth set p = 1/2. To describe strong localization, the Anderson parameter is 
taken to be much smaller than one (V/W <C 1). This corresponds physically to a strongly 
disordered sample where the width of the bands (~ 2V) centered at energies ±W is much 
smaller than their energy difference to the Fermi level. 

The effective hopping matrix element can be computed using a locator expansion [TE| 



|19| . The overlap amplitude (Green's function) between the initial and final sites is given by 

<<£iG(£)ivi> + >=£n/^ (2) 

r i r % - e ir 

where |$) represents the state with a localized electron at the initial site, and |^/ + ) the 
state with a localized electron at the final site; Ef is the Fermi energy which will be set 
to zero, and A is the magnetic vector potential. In principle, the sum is over all paths F 
between the initial and final sites (including back-scattering). However, for V/W <C 1, only 
the shortest (forward scattering) paths need to be included. (For a more detailed discussion 
and justification on this point, using an analogy with high temperature expansions in the 



Ising model, see [fL9| ) P^l .) Neglecting back-scattering, we obtain for paths of length t 

t directed 

w 



(i\G{E)\f) = (L) J( t ), J(t)= £ ]Jvi Tl e iA (3) 

r' i r , 



The sum is now restricted to directed paths T', and T]i = sign (ei) = ±1. The interfer- 
ence information is captured in the function J(t), while the factor (V/W) 1 is the leading 
contribution to the expected exponential decay of the localized wavefunction. 

The transfer matrix approach provides an efficient numerical algorithm for computing 
J{t). As described in reference ||22|| , this method allows summing over the exponentially 



large number of paths in polynomial time (typically ~ t D for D dimensions). The results 
of extensive analytical and numerical studies (mostly in D = 2) based on this method are 
discussed in reference Briefly, the probability distribution for J(t) is quite broad. Its 
logarithm is a universal function with a mean is proportional to t, and variance growing as t 2uJ , 
with u depending on the dimension D. Since the mean and variance of the (log-) distribution 
are independent, two parameters are necessary to describe the tunneling probability. High 



moments of the distribution are however nonuniversal, and dominated by exceptionally good 
samples [ |23| . 

In discussing the change in the tunneling probability in a magnetic field B, NSS in- 
troduced the important concept of the effective 'cigar '-shaped area through which the 
field penetrates. Naively, typical directed paths execute a random walk in the transverse 
direction, so that a path of length t wanders away a distance of the order of t 1 ^ 2 . As shown 
in Fig. || the area presented to a magnetic field perpendicular to such paths thus grows as 
A± oc t x t 1 / 2 oc t 3 / 2 , and the MC is expected to be a function of the flux Bt 3 ! 2 . This is 
indeed the case in the absence of randomness, where the exact response of the sum over 
directed paths in D = 2 |13| initially decreases as B 2 t 3 . The above argument does not work 
in the presence of randomness, where typical paths have super-diffusive transverse fluctu- 
ations that grow as & with ( > 1/2 p2[ . However, the scaling functions are not simply 
modified to depend on Bt l+< *. In the presence of spin-orbit (SO) scattering, the behavior 
is qualitatively similar to the pure case: there is a positive MC, initially scaling as B 2 t 3 , 
which saturates at a finite (t independent) value. By contrast, in the absence of SO, the 
(positive) MC grows unbounded with t. This is because the effect of the magnetic field is a 
(nonuniversal) increase in the localization length £, initially scaling as B 1 ! 2 . The appropriate 
scaling variable in this case is Bt 2 , although numerically one finds a small pre- asymptotic 
regime with Bt 3 ! 2 scaling. There is currently no satisfactory explanation of the crossovers 
in the absence of SO. 



The replica arguments JT9[ suggest that the same asymptotic behavior for the MC should 
be observed in D = 3, as long as the magnetic field is perpendicular to the hopping direction. 
However, it is also possible to consider fields parallel to the hop. The transverse area 
presented to the magnetic field by typical diffusive paths (see Fig. |2|) now grows as Ah oc 
t 1 / 2 x t 1 / 2 , suggesting But as the appropriate scaling argument. This simple argument was 
first presented in ref. [f2(J, along with preliminary numerical support. The anisotropic field 
dependence was verified recently by Lin and Nori P7[ in the IPA approach. In the next 



section we present detailed numerical results pertaining to the anisotropy of MC in D = 3. 
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III. NUMERICAL RESULTS FOR A SINGLE HOP 



As the distribution of the tunneling amplitude J(t) in eq.(^|) is broad, care in averaging 
is quite important. We typically averaged the logarithm of the transition probability (log- 
conductance) over 2000 realizations of randomness. The transfer matrix method allows us 
to examine systems of size t = 600 in the wedge geometry. Furthermore, after studying 
the dependence of the computed amplitudes on the lateral dimension, we also used a bar 
geometry with dimensions 1500 x 200 x 200. This is reasonable if the important paths 
have transverse fluctuations smaller than 200, which was found to be the case. As discussed 
before, our main focus is on the different responses for fields parallel and perpendicular to the 
hop direction. We discuss separately the MC with and without SO scattering. The results 
in the presence of SO are easily interpreted and offer no fundamental surprises. However, 
most of our numerical results in the absence of SO pertain to a pre-asymptotic regime for 
which we have no satisfactory theoretical understanding, but which are most probably of 
experimental importance. 



A. MC without Spin-Orbit Scattering 

Fig. |3| shows the MC and its fluctuations for a magnetic field parallel to the hopping di- 
rection. For the largest values of B, it is clear that ln(| J(t, B)\/\J(t, 0)|) 2 grows linearly with 
the length t of the hop. This is indicative of an exponential correction to the conductance 
due to an increased localization length in the magnetic field. It is only after about t = 400 
that reasonable linearity is achieved, so rather large systems must be examined to study the 
true asymptotic limit. A similar positive MC (and increased localization length) behavior 
is also observed for the perpendicular field orientation. Concurrently, there is a reduction 
in the magnitude of the fluctuations in the tunneling probability (inset Fig. 0), and there 
appear to be strong correlations between changes in the average of the log-conductance and 
its fluctuations. This is also the case in D = 2, where a replica argument suggests that 
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these two quantities (indeed the whole probability distribution) are governed by a single 



parameter 19,22 



Most of the data in Fig. [5] pertains to a pre-asymptotic regime, before the change in 
localization length becomes apparent. Since the length of the hop in most experiments is 
only of the order of 20 to 50 £, it is useful to explore this regime carefully. In Fig. [| we present 
an attempt to collapse the numerical data in this regime for different values of B and t. The 
collapse for the parallel field orientation is demonstrated in Fig. |]a; the maximum hopping 
length in this graph is t — 600, while fields go up to 0.1 flux quanta per plaquette; all in the 
pre-asymptotic regime. The choice of the scaling variable B\\t, is consistent with the flux 
through a section of the NSS 'cigar' perpendicular to the hop direction (A\\ oc \ft x y/t). 
Two regimes are apparent in Fig. |]a: (i) for the lowest fields (5 x 1O~ 5 O per plaquette) and 
sizes (t = 10 — 100) there is a linear increase with the variable Bt. (ii) For intermediate 
fields and hop sizes, when approximately one flux quantum penetrates the NSS 'cigar,' there 
is a non-trivial apparent exponent. The behavior in these regimes is summarized by 

/inffiMV/ 1 -^ f ° ri¥ - 1 f4) 
\ n \J(t,0)\ / \( aiJ B||fr for B ]{ t >l' 

where e*i = 0.38 ± 0.02. 

The corresponding collapse for fields in the perpendicular orientation is presented in Fig. 

||b. In this case the appropriate scaling variable is Bt 3 ^ 2 , again consistent with the flux 

through the NSS 'cigar.' Once more, two different regimes are identified, with 

\ 11 \J(t,0)\ / \ faBJ?/ 2 ) * f or B± t 3 / 2 > 1 ' 
and «2 = 0.25 ± 0.02. We again emphasize that the second regime above is still pre- 
asymptotic. For larger systems (200 x 200 x 1500) the log- conductance crosses over to a 
regime where presumably the relevant scaling variable is Bt 2 [[19] . The latter scaling suggests 



that the magnetic length is the relevant length scale [2^] . We were not able to clearly access 
this regime as cumbersomely large systems must be simulated. 
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It is interesting to compute, on the basis of the above results, the anisotropy in conduc- 
tance of a single critical resistor. We shall define an anisotropy parameter, 

= (ln\J(t,B ± = B)\ 2 -ln\J(t,0)\ 2 ) 
(in \J{t,B\ { = B)\2-ln\J(t,Q)\ 2 Y 
thus making contact with the original experimental definition of ref. 0. Depending on the 
strength of the magnetic field, this anisotropy shows different scaling forms: For the smallest 
fields, such that Bt 3/2 < 1, 

H 1.5Bt 

In this range the anisotropy is field independent, but changes with temperature since t = 
£(T /T) 1//4 . There is an intermediate regime where Bt 3 / 2 > 1 while Bt < 1, and 

P oc b^-H 1 ' 5 "'- 1 , 

which depends on both B and t. Finally, for Bt > 1, 

P oc B a2 - a H 1 - 5a2 ' ai , 

which, using the numerically estimated values, is approximately independent of t, and has a 
weak field dependence as P oc 1/B - 13 . Thus anisotropy is reduced when the field increases 
as shown in Fig. |5|. Qualitatively similar behavior is observed in InO samples for sufficiently 
high fields in ref. ||25|| . 

B. MC with Spin-Orbit Scattering 

When spin-orbit active impurities (doping with heavy elements) are taken into account, 
the NSS model must be generalized to include scattering of the spins. The tight binding 
Hamiltonian is now modified to 

% = , ^2 e i a l,cr a i,o- + X! Vij,<r<r ,a i,(r a j,<r' > (6) 

i,(J <ij>,aa' 

where a indicates the electron spin. The constant nearest-neighbor hopping elements V in 
eq.(|lD are no longer diagonal in spin space. Instead, each is multiplied by Uy, a randomly 
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chosen SU(2) matrix which describes the spin rotation due to strong SO scatterers on each 
bond |T6| , |2"8]j . Eq.(Q) for the overlap of wave-functions at the two end-points must now 
include the initial and final spins, and has the locator expansion 

{ra\G{E)\fa') = Y.\{^ L . (7) 

Each bond along the path contributes a random spin rotation U, and a phase factor from 
the magnetic vector potential A, resulting in 

A = (ta\G(0)\fa') = W(V/Wyj(t); 

with J {t) = Y.Il^ iA U. (8) 

r' i T i 

After averaging over the initial spin, and summing over the final spin, the tunneling proba- 
bility is 

T = Uv(A^A) = W 2 (V/W) 2t I(t); 

with I{t) = hi{J^J). (9) 

For a three dimensional lattice we studied numerically the statistical properties of I{t). 
Using a transfer matrix we evolve paths of length t = 600 in the wedge geometry, averaging 
over 2000 realizations of randomness. As in the previous section, we also used the bar 
geometry to compute I(t) for systems of size 200 x 200 x 1500. We checked the bar-geometry 
results for crossover effects (because of the smaller lateral sizes) and confirmed that their 
behavior is effectively three dimensional. The numerical results are shown in Fig. || and 
have the same qualitative features as in two dimensions: Unlike in the absence of SO, there 
is no linear increase of A In I so (B,t) with system size t, and the MC saturates to a field 
dependent value for sufficiently long hops. The scale of fluctuations in (the logarithm of) 
I(t) is not significantly modified by the B field. As in D = 2, turning on the SO scattering 
from zero, thus switching from an orthogonal to a symplectic Hamiltonian, is accompanied 
by an increase in the zero field conductance, and a concomitant reduction in conductance 
fluctuations. 
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Fig. [7| demonstrates the collapse of the MC data for fields both parallel and perpendicular 
to the hopping direction. For fields parallel to the hopping direction the appropriate scaling 
parameter is again But, corresponding to the flux through the area Au oc (v^) 2 - The scaling 
function has the form 

f 0.4B 2 t 2 if But < 1 
(ln/(t, J B i |)-ln/(t,0))= . (10) 

[C^0.2 if£||i>l 

For fields in the perpendicular orientation, the appropriate area is A± oc t 3 ! 2 (see Fig. [|), 
leading to a scaling function 

(O.IBH 3 tfB ± 2 <l 
(\nI(t,B ± )-\nI(t,0)} = { - - . (11) 

{ C w 0.2 if B ± t 3 / 2 > 1 

The saturation value, when one flux quantum threads the appropriate area, is roughly the 
same in the two cases. 

The anisotropy parameter in the presence of SO is 

f 0.25t if Bt 3 ' 2 < 1 

U if Bt > 1 

with a small crossover region between Bt 3 / 2 ~ 1 and Bt ~ 1. We thus obtain a hopping size 
dependent anisotropy (which might show up as a temperature dependent anisotropy) for low 
fields. For higher fields, anisotropy in the presence of SO scattering disappears as the field 
is increased beyond a flux quantum through the smaller of the typical areas found above. 
In order to properly compare with experiments, one must average over different magnetic 
fields orientations as discussed in the next section. 

The MC in the presence of SO can be explained by a replica analysis. After averaging 
over the impurity potential e« = ±W, the moment (I(B,t) n ), is obtained as a sum over n 
paired paths. The pairings either involve paths taken from J and J^; the contributions from 
the magnetic vector potential cancel for such 'neutral' pairs, which do not contribute to the 
MC. It is also possible to construct pairs with both paths taken from either J or J^; such 
'charged' paths are responsible for magnetic response. An interesting feature of averaging 
over strong SO scatterers (the matrices U in eq.fl8|)) is that the charged and neutral paths 
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become completely decoupled and can be treated independently fL9| . Thus the MC with 
SO is (rather fortuitously) calculated exactly by IPA: there is no change in the localization 
length, and only a constant increase in the tunneling amplitude. Recently Lin and Nori 
p7| computed the field dependencies in a scheme that is equivalent to IPA in two and three 
dimensions. (The equivalence to IPA can be readily recognized from the equality of the 
moments of the 'conductance' to those of a Gaussian obtained in ref. |§.) In the small field 
limit Lin and Nori confirm exactly the scaling variables derived by us here, for the two 
field orientations. Another notable feature of the data is the B 2 dependence for small fields, 
as expected in the IPA approach ||. 

IV. AVERAGING OVER MANY HOPS 

So far we focused on the response of a single hop to a magnetic field. However, the 
percolation arguments [|K]] for the Miller-Abrahams network leads to the conclusion 
that critical resistors dominate only over a correlation length I. Starting from the properties 



of a percolation cluster near the threshold it is concluded |2T[ that a single hop is responsible 
for the overall conductivity only for length scales up to 



where v is the exponent for the divergence of the correlation length close to the percolation 
threshold. The macroscopic system is then built by superposing many blocks of length I. 
Therefore, in general, many critical hops, in general, contribute to the conductivity of a large 
sample. In eq.(|T2|) the variable T /T can be regarded as a measure of disorder strength; it 
increases when the density of states at the Fermi level, or the localization length, decreases. 
As disorder increases, the volume dominated by the critical NSS resistor gets larger. 

For samples of size L ^> £, many hops contribute to the conductivity, and the overall 
MC must be calculated from their average response. Using the experimentally reported |2j 
data, we can estimate £: For samples of thickness d = lOOA and hopping length t ~ 3 — Ad, 
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one obtains t ~ lOOOA. Although this length is large compared to the sample thickness, 



it is still very small in comparison to the sample planar dimensions (10 x 6 mm for the 
experiments in reference ||). Many hops thus influence the conductance of these samples, 
each with a presumably different orientation with respect to the magnetic field. In the 
following paragraphs we shall perform an average over the directions of the NSS 'cigar'- 
shaped' region with respect to the magnetic field. We shall assume that the electric field is 
sufficiently weak so that the relative orientations of the hops are still randomly distributed 
in space p0| , p9| . We expect that for thick samples (as compared to the hopping length 
t) the averaging over all possible directions of the hop removes the anisotropy in MC, as 
seen in experiments. On the other hand, for thin samples such that d < t, the restricted 
averaging over effectively two-dimensional hops should lead to significant field anisotropy. 
Such thickness dependent anisotropy is indeed observed in the experiments of Faran and 
Ovadyahu [@], and those of Orlov and Savchenkov 0,0] ■ 

We shall assume that the appropriate 'cigar' shaped region for calculating magnetic 
interference phenomena is an ellipsoid of revolution as depicted in Fig. |8]. The major and 
minor axes of the ellipsoid are denoted by b and a respectively. Following NSS, 2b = t is the 
length of the dominant hop, while a = y/t£ is a typical diffusive distance in the transverse 
direction. Consider a magnetic field B at an angle 9t with respect to the major axis. We 
take the relevant magnetic flux to be that which penetrates the projection of the ellipsoid 
onto a plane perpendicular to the B field, as indicated in Fig. |8|. This projection is an ellipse 



of minor axis a, and with a major axis of length c = y a 2 cos 2 9 t + b 2 sin 9 t . In previous 
sections we demonstrated that in the weaker field regimes A In | J(B,t)\ 2 oc (|$| 7 ), where $ 
is the appropriate flux. In the presence of SO, 7 = 2 as justified by IPA, while 7 = 1 in the 
absence of SO from the numerical results. In the following, we shall focus on 7 = 2 with 
SO, which has a better justification, and for which it is easier to compute the average, 




($ 2 ) = ((nacBf) 




(13) 
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In performing the average over hop orientations 9t, we distinguish between the following 
regimes: 

(i) The behavior of samples of thickness d 3> t is effectively three dimensional. All orienta- 
tions 9 t are equally likely (in low voltage bias), and (cos 2 9 t ) = 1/3. The averaged response 
is isotropic and depends upon 

($ 2 > = ^V + a 2 ). (14) 

The characteristic B 2 dependence is a signature of the IPA, and is also observed experimen- 
tally for small fields. In fact, the flux through a typical interference region in ref. || is of 
the order of 0.5/ic/e which corresponds to the regime Bt 3 ^ 2 < 1. 

(ii) For a <C d ~ t, the range of orientations of the cigar is limited by the finite thickness, 

leading to MC anisotropy. Consider a magnetic field at an angle a with respect to the plane 

of the sample. We shall indicate the orientation of the major axis of the ellipsoid by a polar 

angle 9 (with respect to the normal to the sample), and an azimuthal angle (p. Because of 

the finite thickness, the range of cos# is limited to the interval [—d/2b, +d/2b], and allowed 

angles in this range are weighted by 

( 2b/d-4b 2 \cos9\/d 2 if |cos#| < d/2b 
p{cos9) = { . (15) 

[0 if\cos9\>d/2b 

The relative angle between the field and the major axis of the ellipsoid is obtained from 

cos 9t = cos a sin 9 cos <fi + sin a cos 9. (16) 

Since all azimuthal angles are possible (cos0) = 0, while (cos 2 0) = 1/2; and from eq.(|15D 
(cos 2 9) = d 2 /8b 2 . Thus, we finally arrive at 

2 . \ d 2 cos 2 a / 3d 2 \ 

cos ')-m + — X -W ' (l7) 



which describes the MC anisotropy, when substituted in eq. (|i~3D . 

(iii) Another limit that is easily accessible is for d ~ a. Now all the ellipsoids lie in the 
plane, and cos# = 0, leading to (cos 2 9 t ) = cos 2 a/2. From eq. (fT3|) we then obtain 
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($2) = [ 26 2 _ (tf _ fl2 ) cog 2 a ] (lg) 

This particular limit is plotted in Fig. || and can be compared to experimental data for 
the angular dependence of MC (at fixed temperature and field strength). The parameters 
a = and 2b = t in this figure are chosen to correspond with those reported in ref. ||. 
The general form of the curve agrees qualitatively with the experimental data. Another 
pertinent comment concerns the experimental data of Laiko et al J7| : they compare results 
for fields parallel and perpendicular to the current direction (while the field is in the plane 
of the sample). On averaging over hop directions, the MC for these configurations should 
be identical, in agreement with experiments, except for very disordered samples where i is 
so large that averaging is not appropriate. 

(iv) Finally for d < a, the NSS cigar is flattened into a pancake and is no longer ellipsoidal. 
The behavior is generally two dimensional, with $ = nab x Bsina. This formula breaks 
down only at very small angles such that a < d/b, for which (<3> 2 ) « (Birda) 2 /2. 

Clearly, the general tendency is that as temperature is reduced (hence 2b = t is increased) , 
the ratio d/t gets smaller. As indicated by the above sequence, this leads to more and more 
pronounced effects in the MC anisotropy. This is indeed consistent with the experimental 
observations. In principle, the temperature dependence of this increase could be measured 
experimentally and compared to the above theoretical formulas. 



V. DISCUSSION 

In this paper we studied the NSS model for quantum interference effects of a tunneling 
electron in D = 3. The three dimensional geometry allows us to consider the relative 
orientations of the hop and the magnetic field. The effect of both spin-orbit active and 
inactive impurities were taken into account. The results indicate that, in the absence of SO 
impurities, there is positive MC due to an increase in the localization length for fields both 
parallel and perpendicular to the hop direction. Furthermore, the MC data for different 
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fields collapses onto universal curves using the scaling variables B±t 3 / 2 , and But for the 
two orientations. This implies that, at least in the low field regime, crossover effects are 
controlled by the flux penetrating an NSS 'cigar' whose transverse size grows diffusively (as 
Vi). 

In the presence of SO impurities, the numerical results again indicate a positive MC, but 
as in two dimensions, no change in the localization length is observed. No reduction in the 
fluctuations is obtained in agreement with the replica arguments ||19|| . The results in this case 
are the same as those obtained from an independent path approximation. The MC grows 
initially as B 2 for both parallel and perpendicular orientations as expected. Its anisotropy 
disappears for large enough fields as the MC saturates to roughly the same constant in both 
directions. This could be checked experimentally in the thin film limit (d < t) for samples 
doped with heavy elements. 

The most spectacular manifestation of these results is the possibility of observing 
bulk MC anisotropy, when the sample is small enough (or the characteristic length I = 
^(To/T)^ +1 ' ) ^ D+1 ^ is large enough) that only a single hop (or just a few) dominates the 
macroscopic conductance. The sample will then exhibit a 'finger-print' in its MC anisotropy, 
and the random orientation of the critical hop can be determined by bulk anisotropy mea- 
surements. (It is possible that this phenomenon explains why the experiments of Laiko et 
al JT] go from isotropic to anisotropic behavior as disorder is increased fl5C[l .) Large values 
for £ may be achieved by either choosing samples with lower density of states at the Fermi 
level, or larger Other manifestations of the disorder length scale (including bulk MC 



anisotropy) may occur under strong voltage bias |[29|| , if the electric field reorients the critical 
hops; constraining the average over hop orientations to within a cone. 

However, the conductance of most samples with L 3> i is governed by many hops. Our 
estimate of I based on published experimental data || indicates that it is much smaller than 
typical sample planar dimensions. Thus, averaging over many hop orientations is inevitable, 
washing out the predicted anisotropy finger-print of a single hop. Nevertheless, anisotropic 
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behavior is still expected for sufficiently thin samples (or at low temperatures). This is 
because when the length of the cigar becomes comparable or larger than sample thickness 
(t > d), the hops are forced to lie mostly parallel to the sample plane. This restriction on hop 
orientations then leads to an MC anisotropy that becomes more and more pronounced upon 
lowering temperature. Appearance of such anisotropy in thin films has already been observed 
in insulating InO samples |25|||. In principle, the variations of anisotropy with temperature 
can be measured accurately and compared to the formulas derived in the previous section. 
Measurements of anisotropy can thus provide an additional experimental tool for tests of 
quantum interference models. 
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FIGURES 

FIG. 1. The NSS model on a three dimensional diagonal lattice. Impurities are located on the 
sites. 

FIG. 2. The figure depicts how effective areas that arise from scaling can be derived from 
random walk arguments. 

FIG. 3. Log-conductance as a function of system size t for fields parallel to the hop direction. 
The change in the slope with the field indicates an exponential correction to the conductance 
(change in the localization length). A straight line is drawn as a guide to the eye. The inset shows 
a reduction of fluctuations with the field. The power law dependence on the hopping length is also 
indicated. 

FIG. 4. The figure shows the collapse of the log-conductance data with the appropriate scaling 
variables (a) Bt in the case the field is parallel to the hop and (b) Bt 3 ^ 2 for fields perpendicular to 
the hop direction. The scaling variable tells about the relevant area threaded by the field. 

FIG. 5. Field dependencies for parallel and perpendicular field directions when Bt 3 / 2 > 1 and 
for a single critical resistor. 

FIG. 6. Magnetoconductance in the presence of spin-orbit scattering. The relevant 
log-conductance is no longer linear in t for large hopping lengths, indicative of no changes in 
the localization length due to the field. 

FIG. 7. Collapse of the MC data in the presence of SO scattering. The MC is governed by the 
areas a ) (^/H 1 / 2 ) 2 for fields parallel to a single critical hop and b) £}/ 2 t 3 / 2 for fields perpendicular 
to the hop direction. 

FIG. 8. Angle convention used to perform averages over critical hop directions. The magnetic 
field is along the z axis making an angle 9t with the principal axis of the ellipsoid containing 
dominant paths. The relevant magnetic flux penetrates through the largest cross section of the 
ellipsoid perpendicular to the magnetic field i.e. an ellipse of minor axis a and major axis c. 

23 



FIG. 9. Theoretical curve computed from eq.(18). The experimental data correponds to 
£ = 85A, t = 280A and sample thickness d = 250A Notice how the curve fall towards a = 
(parallel to sample) , closely resembling the experimental data of reference [2] . 
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Coherent hop in the Miller-Abrahams network 
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